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LONG-TERM  GOALS 

The  long-term  scientific  goal  of  my  research  is  to  better  understand  the  distribution  of  phytoplankton 
in  the  world's  oceans  through  remote  sensing  their  influence  on  the  optical  properties  of  the  water.  An 
associated  goal  is  the  understanding  of  the  absorption  and  backscattering  properties  of  marine  particles 
in  terms  of  the  distributions  of  their  size,  shape,  and  composition. 

OBJECTIVES 

The  inherent  optical  properties  (lOPs)  of  marine  particles  are  most-often  modeled  as  homogeneous 
spheres  using  Mie  Theory.  Although  this  approach  has  been  fruitful,  the  next  logical  step  in  modeling 
marine  particles  is  to  abandon  the  normally-employed  spherical  approximation  and  use  more  realistic 
approximations  to  their  shape.  The  advent  of  computer  codes  capable  of  handling  more  complex 
shapes,  and  the  increased  computational  speeds  now  available,  suggest  that  particle  modeling 
employing  simple  non-spherical  shapes,  e.g.,  disks,  rods,  etc.,  could  become  routine.  For  example, 
Gordon  and  Du  (2001)  used  a  two-disk  model  to  try  to  reproduce  the  backscattering  by  coccoliths 
detached  from  E.  huxleyi  in  the  blue-green  region  of  the  spectrum.  They  found  that  the  gross 
morphology  of  the  particle  (e.g.,  modeling  the  coccolith  as  parallel  disks  as  opposed  to  a  single  disk  of 
the  same  diameter  and  volume)  was  paramount  in  determining  the  spectral  variation  of  backscattering. 
Later  Gordon  et  al.  (2009)  showed  that  the  resulting  spectral  variation  of  the  backscattering  cross 
section  for  the  two-disk  model  agreed  with  experiment,  and  that  additional  fine  structure  evident  in 
SEM  images  of  coccoliths  had  only  a  minor  influence  on  the  backscattering.  In  another  study,  Clavano 
et  al.  (2007)  looked  at  the  lOP  of  spheroids  (ellipses  of  revolution)  as  a  function  of  their  size  and 
aspect  ratio  (larger  major  axis  divided  by  smaller  major  axis)  using  refractive  indices  in  the  range 
found  for  marine  particles.  Their  results  suggest  that,  when  compared  to  equal  volume  spheres,  shape 
can  be  a  significant  factor  in  the  lOPs,  especially  backscattering. 

My  present  work  is  a  continuation  of  my  research  into  computation  of  the  lOPs  of  non-spherical 
objects.  Of  particular  interest  to  me  is  the  scattering  and  absorption  of  particles  showing  significant 
deviation  from  spheres,  e.g.,  cylinders  with  large  aspect  ratios.  This  gross  morphology  resembles  long 
chain  phytoplankton.  How  do  their  scattering  and  absorption  properties  depend  on  the  length,  diameter 
and  index  of  refraction  of  the  cylinders,  or  more  significantly,  their  aspect  ratio  (AR  =  length  -h 
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diameter)?  Do  their  lOPs  reaeh  a  limit  in  whieh  they  beeome  simply  proportional  to  the  eylinder’s 
length,  so  their  length-normalized  lOPs  beeome  independent  of  aspeet  ratio?  How  do  they  depend  on 
the  distribution  of  absorbing  pigments  within  the  eylinder?  How  do  they  depend  on  orientation?  How 
do  they  eompare  to  the  seattering  and  absorption  by  spheres  eontaining  the  same  volume  of  material 
and  quantity  of  absorbing  pigment?  Is  the  spherieal  approximation  effeetive  in  interpreting  their 
seattering  and  absorption  in  terms  of  an  equivalent  refraetive  index?  My  goal  is  to  develop  an 
understanding  at  a  level  that  would  allow  the  ineorporation  of  shape  information  into  lOP  models  and 
in  the  analysis  of  lOP  data. 

APPROACH 

There  are  now  several  eomputer  eodes  that  allow  aeeurate  eomputation  of  the  seattering  and  absorption 
properties  of  individual  partieles  with  smooth  nonspherieal  shapes.  I  have  used  the  diserete-dipole 
approximation  (DDA,  Draine,  1988;  Draine  and  Flatau,  1994)  in  my  study  of  partieles  similar  to  E. 
huxleyi  eoeeoliths.  I  am  now  using  a  later  version  of  the  eode  (DDSCAT  7.0)  to  eompute 
baekseattering  by  randomly  oriented  finite  ey finders  with  length  up  to  25  pm  and  diameters  up  to  1.5 
pm.  At  this  size,  an  equal  volume  sphere  has  a  diameter  of  about  4  pm. 

WORK  COMPLETED 

I  have  eompleted  a  study  of  the  seattering  and  baekseattering  of  non-absorbing  randomly-oriented 
eylinders  (diameter  of  0.5- 1.5  pm)  as  a  funetion  of  the  length  and  the  index  of  refraetion  (1.02  to  1.20). 
The  goal  was  to  understand  at  what  aspeet  ratio  the  baekseattering  beeomes  essentially  proportional  to 
the  length.  This  was  mostly  earried  out  using  the  DDSCAT  6.1  eode.  For  the  study  of  absorbing 
partieles,  I  have  modified  the  new  DDSCAT  7.0  eode  to  inelude  eylinders  with  internal  strueture,  i.e.,  a 
eoated  eylinder.  This  is  ehosen  to  represent  a  long  ehain  of  phytoplankton,  e.g.,  diatoms,  in  whieh  the 
absorbing  pigment,  if  any,  is  paekaged  within  the  internal  eylinder. 

I  have  also  eonsidered  absorbing  eylinders  ranging  in  diameter  from  0.5  to  1.5  pm  and  lengths  from 
0.5  to  15  pm.  Two  sets  of  eylinders  I  examined  were  homogeneous  and  had  refraetive  indiees 
1.05-0.002/  and  1.05-0.010/.  Two  other  sets  were  eoated  eylinders  in  whieh  the  inner  eylinder  had  a 
diameter  equal  to  half  the  eylinder  diameter  and  indiees  of  1.05-0.008/  and  1.05-0.040/.  The  outer 
eladding  had  an  index  of  1.05  (no  absorption).  The  eylinder  with  homogeneous  index  of  1.05-0.010/ 
and  the  eoated  eylinder  with  eore  index  1.05-0.040/  have  exactly  the  same  quantity  (mass  or  number 
of  molecules)  of  absorbing  pigment. 

The  quantities  that  I  computed  are  the  orientationally-averaged  extinction  cross  section,  (Tc,  the 
scattering  and  baekseattering  cross  sections,  o*  and  (Jhb,  the  absorption  cross  section,  (Ja,  and  the 
scattering  phase  function  (differential  scattering  cross  section  divided  byo/,).  The  extinction  and 
absorption  efficiencies,  Qc  and  Qa,  are  the  associated  cross  sections  divided  by  the  orientationally 
averaged  projected  area  of  the  particle  (the  area  of  the  particle’s  shadow  formed  by  the  incident  beam, 
LD/2). 
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RESULTS 


The  main  finding  of  the  study  is  that  for  AR  >  ^  3-5,  Qc  and  Qa,  beeome  independent  of  AR  and,  in 
faet,  elose  to  those  computed  for  infinite  cylinders.  Similar  results  are  obtained  for  the  scattering  phase 
function  and  the  backscattering  probability  {cThtJcTh)- 

Figure  1  provides  the  computed  values  of  Qc  and  Qa  for  the  all  aspect  ratios  examined  with  absorbing 
cylinders  {mi  >  0)  and  mr=  \  .05.  The  notation  in  the  legend  is  provided  in  the  figure  caption.  The 
upper  two  panels  are  for  all  aspect  ratios  (1/3  -  30)  and  the  lower  panels  for  AR  >  3.  The  lines 
correspond  to  exact  computations  for  infinite  cylinders  with  the  same  diameter  computed  with  the 
IPHASE  code  (Evans,  1988;  Eournier  and  Evans,  1996).  Eor  the  case  with  weaker  absorption,  the 
results  clearly  show  that  the  efficiencies  closely  follow  those  for  infinite  cylinders  as  long  as  AR  >  3, 
for  both  the  coated  and  the  homogeneous  cylinders.  Eor  the  case  with  stronger  absorption,  the  so- 
called  “package  effect,”  the  decrease  in  absorption  when  the  absorbing  molecules  are  not  uniformly 
distributed  within  the  particle,  results  in  lower  absorption  efficiency. 

This  packaging  effect  displayed  more  clearly  in  Eigure  2,  for  which  we  provide  the  ratio  of  the 
absorption  efficiencies  Qa  for  a  coated  cylinder  with  indices  =1.05-0.040/  and  nioutside 
1.05-0.000/  (ga(Packaged))  to  that  of  a  homogeneous  cylinder  nimside  =  ^outside  =  1.05-0.010/ 
(ga(Homo)).  Note  that  both  cylinders  contain  the  same  number  of  absorbing  molecules.  The  figure 
shows  that  the  effect  of  the  absorbing  pigment  packaging  is  greatest  in  the  blue  region  of  the  spectrum 
and  for  larger-diameter  cylinders.  The  maximum  decrease  in  Qa  due  to  the  packaging  is  about  25%  in 
this  case.  Although  the  symbols  do  not  differentiate  between  cylinder  lengths,  for  a  given  diameter  the 
packaging  effect  is  largest  in  the  shortest  cylinder  and  depends  very  little  on  the  length  once  AR 
exceeds  unity.  In  the  case  with  less  overall  absorption,  i.e.,  minside  =1.05-0.008/  and  nioutside  = 
1.05-0.000/  compared  to  that  of  a  homogeneous  cylinder  =  nioutside  =  1.05-0.002/,  similar 
results  are  obtained  (not  shown);  however,  the  maximum  decrease  in  Qa  due  to  the  packaging  is  only 
about  10%. 

The  extinction  efficiencies  of  finite  cylinders  with  larger  values  of  p  are  compared  with  those  of 
infinite  cylinders  in  Eigure  3.  In  this  case  the  refractive  index  is  1.20,  and  since  there  is  no  absorption, 
Qb  =  Qc  ■  The  figures  clearly  show  that  iox  AR  >  3,  the  scattering  efficiency  is  again  close  to  that  of  an 
infinite  cylinder,  with  all  cases  except  two  differing  from  an  infinite  cylinder  by  less  than  ±10%  (RMS 
difference  ~  5%). 

We  have  seen  that  the  extinction  and  absorption  efficiencies  of  micrometer-sized  cylindrical  particles 
depend  little  on  the  aspect  ratios  as  long  as  AR  >3.  Is  this  the  case  for  the  phase  function  and 
backscattering  probability?  Eigure  4  provides  an  example  of  the  weak  dependence  of  the  scattering 
phase  function  (and  degree  of  linear  polarization)  on  aspect  ratio.  Virtually  the  only  differences  in  the 
phase  function  between  AR  =  20  and  AR  =  3  are  the  enhanced  scattering  near  zero  degrees  and  the 
deeper  minima  near  30°,  50°,  80°,  and  130°  iox  AR  =  20.  This  weak  dependence  on  aspect  ratios  as 
long  as  AR  >  5  is  also  displayed  by  the  backscattering  probability  as  shown  in  Eigure  5  for  a  wide 
range  of  refractive  indices  and  particle  diameters. 
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IMPACT/APPLICATIONS 


The  fact  that  (normalized)  lOPs  by  cylinders  with  aspect  ratios  of  3-5  become  independent  of  the 
aspect  ratio  can  be  used  to  simplify  computation  of  such  for  analysis  of  particle  scattering.  For 
example,  computation  of  the  scattering  efficiencies  for  infinite  cylinders  (which  can  be  carried  out  very 
efficiently)  can  be  used  to  estimate  the  same  for  cylinders  with  AR  >  3  with  an  error  <  5-10%.  These 
efficiencies  can  then  be  used  to  compute  the  cross  sections  for  the  cylinder  length  of  interest.  Similarly, 
the  backscattering  probability  for  long  cylinders  with  diameters  -'1-3  times  the  wavelength  (which 
cannot  be  computed  with  most  light  scattering  codes)  can  be  deduced  (within  5-10%)  from  those  of 
finite  cylinders  with^fi?  -  3-5.  These  observations  can  simplify  the  analysis  of  light  scattering  and 
absorption  by  high-aspect  ratio  particles. 

RELATED  PROJECTS 

None. 
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Figure  1.  Extinction  and  absorption  efficiencies  computed  for  randomly  orientated,  homogeneous 
or  coated  cylindrically  shaped  particles,  given  that  their  diameter  and  aspect  ratio  are  known  so  that 
their  orientationally-averaged projected  area  is  DL/2.  Here,  p  =  2a  (my-l)  and  p'=  4a  ffi/,  where 
mr  -imi,  is  the  refractive  index  of  the  particle  relative  to  water,  and  a  =  ttD/A,  with  D  the  cylinder’s 
(outer)  diameter  and  X  the  wavelength  of  light  in  the  water.  Solid  lines  are  the  exact  computations 
for  randomly  oriented,  infinite  cylinders.  The  notation  “m=1.05-040i_1.05-000i”  indicates  that  the 
refractive  index  of  the  core  is  1.05  -  0.040i,  and  the  refractive  index  of  the  coating  is  1.05  -  O.OOOi, 
etc.  In  the  case  of  coated  cylinders,  p  and  p'  are  computed  using  mi  of  the  associated  homogeneous 
particle.  Top:  all  aspect  ratios  (1/3  -  30).  Bottom:  all  aspect  ratios  >3.  The  graphs  clearly  show  that 
the  efficiencies  closely  follow  those  for  infinite  cylinders  as  long  as  AR  >3,  for  both  the  coated  and 
the  homogeneous  cylinders.  For  the  case  with  stronger  absorption,  the  so-called  “package  effect,  ” 
the  decrease  in  absorption  when  the  absorbing  molecules  are  not  uniformly  distributed  within  the 

particle,  results  in  lower  absorption  efficiency. 
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Figure  2.  This  figure  provides  the  ratio  of  absorption  efficiencies  (coated  to  homogeneous)  of  the 
more  strongly  absorbing  cylinders  as  a  function  of  wavelength.  The  diameter  of  the  cylinder  (in  pm) 
is  specified  in  the  legend.  For  each  diameter,  the  symbols  refer  to  cylinder  lengths  ranging  from  0.5 
to  15  pm.  The  figure  shows  that  the  effect  of  the  absorbing  pigment  packaging  is  greatest  in  the 
blue  region  of  the  spectrum  and  for  larger-diameter  cylinders.  Although  the  symbols  do  not 
differentiate  between  cylinder  lengths,  for  a  given  diameter  the  packaging  effect  is  smallest  in  the 
shortest  cylinder  and  depends  very  little  on  the  length  once  the  aspect  ratio  (length/diameter) 

exceeds  unity. 
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Figure  3.  Qc  computed  for  non-absorbing  cylinders  (m  =  1.20  -  O.OOOi)  with  diameters  (D)  ranging 
from  0.5  /jm  to  2.0  pm.  Left:  0.25  <AR  <30  (points  colored  in  red  areAR  =  2).  Right:  3  <AR  < 
30.  The  solid  curve  is  the  extinction  efficiency  (for  a  unit  length)  of  randomly-oriented  infinite 
cylinders.  As  in  Figure  1,  p  =  2a(m -1)  with  a=  ttD/A.  The  figures  clearly  show  that  for  AR  >3, 
the  scattering  efficiency  is  again  close  to  that  of  an  infinite  cylinder,  with  all  cases  except  two 
differing  from  an  infinite  cylinder  by  less  than  ±10%  (RMS  difference  ~  5%). 


0  (Deg.)  ©  (Deg.) 


Figure  4.  Orientationally-averaged  scattering  phase  functions  (left)  and  degree  of polarization 
(right)  at  600  nm  (vacuum)  for  homogeneous  cylinders  with  a  diameter  of  1  pm  and  length  of  3  pm 
(DlxL3)  and 20  mm  (DlxL20).  The  refractive  index  is  1.05  -  0.0021  Virtually  the  only  differences 
in  the  phase  function  between  AR  =  20  and  AR  =  3  are  the  enhanced  scattering  near  zero  degrees 

and  the  deeper  minima  near  30  °  50  °  80  ° 
and  130  °for  AR  =  20. 
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Figure  5.  Examples  of  the  variation  of  the  backscattering  probability  with  aspect  ratio  for  cylinder 
diameters  between  0.5  and  1.5  fjm  and  refractive  indices  ranging  from  1.02  to  1.20.  The  black 
curves  are  for  a  vacuum  wavelength  of 400  nm  and  the  red  curves  for  700  nm.  The  graphs  show 
that  there  is  very  little  variation  in  the  backscattering  probability  with 

aspect  ratio  for  AR  >  5. 
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